Measurements of the density-dependent many-body electron mass 
in 2D GaAs/AlGaAs Heterostructures 
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We determine the density-dependent electron mass, m*, in two-dimensional (2D) electron sys- 
tems of GaAs/AlGaAs heterostructures by performing detailed low-temperature Shubnikov deHaas 
measurements. Using very high quality transistors with tunable electron densities we measure m* in 
single, high mobility specimens over a wide range of rs (6 to 0.8). Toward low-densities we observe 
a rapid increase of m* by as much as 40%. For 2 > rs > 0.8 the mass values fall ~ 10% below the 
band mass of GaAs. Numerical calculations are in qualitative agreement with our data but differ 
considerably in detail. 

PACS numbers: 71.18.-(-y, 73.40.-c, 73.43.Qt 



In a crystal the mass of an electron often differs from 
its mass, mo, in free space. The electron mass in a semi- 
conductor can deviate from toq by more than a factor of 
ten. The origin of this effect is interference of the electron 
wavefunction with the periodic array of the ions in the 
solid. Such "single particle" effects are well understood 
and readily calculated. However, there are other factors 
that affect the electron mass. In general, any excitation 
of the solid ~ such as phonons, spin waves, plasmons - can 
impact the dispersion relation of the carrier , IJ, but only if 
such excitations come close to resonance. More intricate 
yet, the electron mass is modified by interactions with 
all neighboring conduction electrons. Naively one would 
think that interactions always enhance the carrier mass 
since they imply pushing against other electrons, making 
them apparently heavier. Yet, curiously, theory tells us 
that the mass can be reduced as weliy. Such interac- 
tions are of complex "many particle" origin and bring us 
to the edge of the theoretical and numerical abilities in 
condensed matter theory. 

The impact of electron-electron (e-e) interaction on 
carrier mass increases on lowering the dimensionality. 
Three dimensional (3D) electron systems are expected 
to show little effect, whereas ID systems are highly af- 
fected. From an experimental point of view, 3D systems 
are most abundant and best characterized, whereas ID 
systems are rare and suffer from many complications such 
as Peierls instabilities and their sensitivity to disorder. 
Two-dimensional electron systems (2DESs) provide an 
excellent compromise for the study of such many-body 
phenomena: 2DESs have been honed to extremely high 
quality 3] and the expected effects are moderately strong. 
In addition, the electron density in a 2DES is in principle 
continuously tunable, allowing the study of such phenom- 
ena within a single specimen over a wide range of densi- 
ties. The strength of e-e interaction is usually described 
by a dimensionless parameter rs, which is defined as the 
ratio of Coulomb interaction to Fermi energy. In 2D, r^ is 



inversely proportional to the square root of density, so a 
variable density translates into a tunable e-e interaction. 

Starting in the late sixties^, there have been many 
theoretical studies '3| of the impact of e-e interaction on 
carrier mass using ever more powerful numerical tools. 
Experimentally, the Silicon Metal-Oxide-Semiconductor- 
Field-Effect- Transistors (MOSFETs) had been the domi- 
nant implementation of a 2DES. Smith and StileS|S| were 
the first to measure m* as a function of in such a de- 
vice. However, subsequent work by Fang et al.0 asserted 
that such mass measurements in MOSFETs were affected 
by several side effects. Yet, the Smith and Stiles 'F] data 
remain the experimental reference point for this rapidly 
progressing area of theoretical investigation. 

Recently, Coleridge et al.0 have performed mass 
measurements on five fixed density GaAs/AlGaAs het- 
erostructures, but only for r^ < 1.7. They observed a 
monotonically increasing m* with increasing carrier den- 
sity. With the recent interest in a transition from an elec- 
tron liquid to an electron solid at very high r^, several 
new studies have taken placed, 0. They concentrate on 
phenomena associated with the anticipated phase tran- 
sition but not with the electron liquid per se. Given the 
theoretical progress in the area of e-e interactions and the 
availability of very high-quality 2DESs, a careful mea- 
surement of the effective 2D electron mass over a wide 
range of density seems to be of considerable importance 
to make contact between theory and experiment. 

Toward this end we have measured m* in a very high 
quality, tunable density, GaAs/AlGaAs heterojunction- 
insulated gate field effect transistor (HIGFET) for 6 > 
Ts > 0.8. We observe a strongly increasing mass towards 
low densities and a mass ^ 10% below the band mass for 
all 2>rs> 0.8. 

Our primary sample, HIGFET- 1, was grown by molec- 
ular beam epitaxy (MBE) onto a (001) GaAs substrate. 
The 2DES resides at the interface of 5/im of GaAs and 
5nm of AlAs, topped with 4/im of Alo.33Gao.67As. The 
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FIG. 1; Shubnikov deHaas (SdH) oscillations for n — 
5.4 X lO^Vcm^ at T=101, 180, 260, 380mK. The lower in- 
set: Fourier spectrum from the lOlmK data vs 1/B. Upper 
inset: all Dingle plots for seven different T ranging from 100 
to 400mK data sets yield similar scattering time TqS. 

latter acts as an insulator separating the channel from 
the gate, which consists of a 25nrn thick heavily doped, 
conducting GaAs ri+ layer. The material was processed 
into a 600/im square mesa using photolithography and 
contacted via Ni-Ge-Au annealed pads. An extra pad 
contacted the gate. Another sample, HIGFET-2, dif- 
fers from HIGFET-1 in a thinner channel material (2/im 
GaAs) and the lack of the thin AlAs layer. Both samples 
have a peak mobility of 1 x \{f cm? /V s. The electron 
density of the 2DEG can be tuned by changing the gate 
voltage of the transistor. This provides a major advan- 
tage as compared to fixed density specimens since the 
strength of e-e interaction can be changed continuously. 

Measurement of the specific heat would be the most 
direct way to determine the mass that includes e-e inter- 
actions. Such experiments are exceedingly challen gin g 
and have not yet been realized with high precision [lOj. 
A mass determined by cyclotron resonance will not ex- 
hibit effects from e-e interactions according to Kohn's 
theorem^J- Instead, we employ the Shubnikov dcHaas 
(SdH) effect and determine the effective mass from the 
reduction of the amplitude of these magnetoresistance 
oscillations with increasing temperature. For SdH mea- 
surements to provide reliable mass data, data collection 
has to be performed within appropriate parameter win- 
dows, and the interpretation of the data has to be con- 
ducted with considerable care and multiple cross checks. 
In the following paragraphs we detail our procedure. 

All measurements were performed in a dilution refrig- 
erator over a temperature range from 100 to 400mK. We 
used conventional low frequency AC lock-in techniques 
with excitation currents ranging from InA to lOOnA, cho- 



sen to avoid sample heating. Data were taken in a single 
cool down in order to have the most consistent quality 
from each sample. The temperature is based on a cal- 
ibrated Ruthenium oxide thermometer mounted on the 
same silver sample holder. Magneto-resistance is negli- 
gible at low fields and the relatively high temperature 
ensures good thermal coupling between sample and ther- 
mometer. At each density, set by the gate voltage, we 
recorded a family of SdH traces at a range of tempera- 
tures. Data collection was limited to the moderate mag- 
netic field region such that the SdH oscillations were well 
developed but before the small spin splitting in GaAs was 
resolved. Therefore, in the region of our measurements, 
each minimum in the oscillations corresponds to a Lan- 
dau level index i. The temperature was kept sufficiently 
high to avoid the quantum Hall regime, in which SdH 
oscillations are becoming non-sinusoidal. 

SdH oscillations are a result of the comb of Landau 
levels sweeping through the Fermi level while the mag- 
netic field is ramped. Hence, the oscillating part of the 
magneto resistivity l^pxx can be written as a Fourier 
series: [T^ITsj . 

^Pxx{e) = Po V'7tftCa,pexp( — ^) cos[27rp(-^ - i)], 
p w*T, hw* 2 

(1) 

with ljj* — eB/m*, and Ca,p being a disorder 
coefficient The T-dependence "fth is given by 

^ p ■ 2TT^kBT/hw* 

sinhip -271^ kBT/hiu*) ^ ' 

It is common practice to maintain only the fundamental 
term and neglect all higher order Fourier components. 
We will later explicitly check this assumption for our 
measurements. We derive the effective electron mass, 
m*, by fitting, on a log(Ai?/T) vs T plot, expression (2) 
to our T-dependent data from each index i, where R is 
resistance. We achieve excellent fits to all SdH data with 
a single mass value for each density in the whole tempera- 
ture range lOOmK < T < AOOmK JA\- Before discussing 
these results we performed several cross-checks to estab- 
lish the reliability of our data reduction. 

In recent literature there arose a concern as to whether 
the 3D SdH formalism applies correctly to 2D cases [T5l|. 
Significant deviations from the original Lifshitz and Kose- 
vich (LK) formula were observed in de Haas-van Alphen 
effects of layered organic conductors'^ as well as in 
2DESs of GaAs|T3. Such deviations can be traced back 
to significant contribution from higher order Fourier com- 
ponents. When higher harmonics are negligible, the 
thermal reduction of the amplitude is well described by 

In order to explore the variability of our mass data due 
to such possible deformations, we have used three differ- 
ent methods to derive m*. i) Data points depicted by 
solid black circles in fig. 2(a) are m* values fitted directly 
to eq.(2) with p=l only. Each point represents the av- 
erage mass from several different Landau level indices, 
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FIG. 2: Density dependence of 2D electron effective mass 
m* (a) m* from three different fitting methods, from a sec- 
ond specimen, and from another groupQ, see legend. Mass 
enhancement due to non-parabolicity shown as dashed line, 
(b) m*s of "fundamental fit" in (a) corrected for nonparabol- 
icity. Several theories are indicated as lines, see legend and 
text. 



and this root-mean-square deviation dominates the er- 
ror bars, ii) Masses depicted as solid squares are derived 
from the same data, but after Fourier filtering them on 
a 1/B plot. A typical Fourier spectrum is shown in the 
inset of fig.l, demonstrating that higher harmonics are 
practically negligible. Accordingly, the filtered and un- 
filtered data differ only slightly from one another, iii) 
We follow a revised version of the LK derivation for 2D 
Fermi liquids under strong e-e interaction and electron- 
impurity scattering; 20]. There, the thermal reduction of 
the first Fourier component of the SdH amplitude follows 
an exponential decay form instead of eq.(2). The masses 
obtained from these fits are shown as crosses in fig. 2(a). 
These m* values are close to those from the Fourier fil- 
tered method. We conclude that our SdH data fall within 
the window over which an LK analysis can be applied and 
that our mass derivation is robust. 

Another possible source of concern are density inho- 
mogeneities. Such inhomogeneities were shown to con- 
siderably affect the determination of the scattering time 



via SdH|2l|; a subject we will address later in this pa- 
per. From the known density gradient across the wafer, 
we determine the density gradient to be < 0.07% across 
our sample. The contribution of density inhomogeneities 
from surface roughness at the interface is estimated to 
be < 0.1% 22). Our analytical calculations and many 
numerical simulations demonstrate that the thermal re- 
duction factor is not affected by such density inhomo- 
geneity and the m* readings at any fixed field will not be 
influenced. 

Data from HIGFET-2 are shown as stars in fig. 2(a). 
Despite the differences between the two samples, their m* 
coincide within the error bars. We also include data from 
four fixed density samples measured earlier by Coleridge 
et al.0 at higher densities. These data are also consistent 
with our results. The m* data from Si-MOSFETs are 
generally - 20% larger than GaAs data.[i|ll3 

The effective mass versus density data of fig. 2(a) fol- 
lows a smooth but non-monotonic curve. With increasing 
Ts at the low density part, there is a strong enhancement 
of m*. At ~ 5, m* exceeds the GaAs band mass by 
^ 40%. In the high density region, m* stays ~ 10% 
below the band mass for 2 > > 0.8. In this regime 
the non-parabolicityj2^ of the GaAs conduction band 
actually further enhances the discrepancy between single 
particle mass and our measured mass. The dashed line 
shows the result of numerical calculation for the band 
mass at the Fermi energy due to 2D confinement plus 
band filling. Results from three different trial wave func- 
tions fall within the thickness of the line. The ratio of 
the measured m* to the band mass, corrected for this 
non-parabolicity, arc shown in fig. 2(b) together with the- 
oretical calculations. For clarity, we limit the data points 
to a subset of the data of fig. 2(a). 

Kwon, Ceperley, and Martin'24'|(KCM) using a varia- 
tional Monte Carlo method, and Asgari et al. [25l| fAfcal) 
using a many-body local field approach have performed 
extensive theoretical calculations of the impact of e-e in- 
teractions on the mass in 2D. The results of these nu- 
merics are plotted in fig. 2(b) as a full solid line and as 
a dashed line, respectively. The KCM theory reproduces 
quite well the average mass values at low , whereas the 
A&al(0) theory seems to depict the overall shape and the 
upturn at high r^. However, both theories assume a zero- 
thickness 2DES. When adjusted for finite thickness, the 
A&al result follows the dotted fine, A&al(f)[^. No such 
adjustment is available for KCM and the coincidence of 
both curves must be considered accidental. In any case, 
these calculations seem to capture some aspects of the 
data, but clearly do not describe the high regime. 

All numerical calculations have been performed for a 
2D system in the absence of a magnetic-field. Since our 
data were collected in small B-fields and mostly in very 
high Landau levels we can regard our data as representing 
this limiting B = case. However, Smith, MacDonald 
and Gumbs 4] (SMG) have performed mass calculations 
in the presence of a B-field based on the random-phase 
approximation (RPA). Their results are shown as a se- 
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quence of thin lines identified by the representative Lan- 
dau index i. We note that i = oo is equivalent to B — 
and can be directly compared to the other calculations, 
showing further the discrepancies between different the- 
ories. These i — oo results also differ considerably from 
our data. At the same time, the SMG calculation shows 
a considerable dependence of the mass value on the index 
which needs to be taken into account. 

The experimental window dictated by temperature, si- 
nusoidal lineshape and available B-field, results in a cor- 
relation between density n and available index i. Lower 
densities require the recording of low indices whereas 
higher densities allow to measure much higher indices. 
While the i = oo SMG results differ considerable from 
our data, it appears reasonable to use their fractional de- 
pendence of the mass on index, m*{i = oo)/m*(?), to 
correct our data for such an index dependence. The re- 
sulting m* data are shown as filled diamonds in fig. 2(b). 
The overall shape of the density dependent mass is not 
much affected although the very low density masses are 
enhanced beyond the previous error bars, since such data 
were taken at relatively low i. On the other hand the high 
density data are almost unaffected. While the correc- 
tion suggested by the SMG work modifies somewhat the 
comparison between theory and experiment, the overall 
conclusions remain intact. 

In addition to m*, the envelope of the SdH oscillations 
provides us with a measure of the quantum scattering 
time Tq at each temperature. Since Tq enters the LK ex- 
pression, eq.(l), a temperature dependence of Tq could 
affect the value of m*. Before examining the data, we 
should stress that a T-dependence of Tq must be con- 
sidered very weak and only on the scale of T/Tp, Tp 
being the Fermi temperature, since scattering by fixed 
imperfections is the only remaining mechanism and it is 



practically T-independent in our temperature and den- 
sity range. Nevertheless, we evaluated Tq employing the 
semi-log Dingle plot of the SdH amplitude normalized 
to po7th(eq.(l)) vs 1/B. The upper inset of fig.l shows 
such a typical Dingle plot for seven different tempera- 
tures. The slope of the data determines Tq, which can be 
affected by density inhomogeneities[2l| . However, here 
we pursue only a possible T-dependence of Tq. From ex- 
tensive modeling we determine that it is unaffected by in- 
homogeneity of as much as 10%. Evaluating many Dingle 
plots, we find that over our T-range, (T) varies by less 
than 1% for all rg < 4, see also ref.Q- At the lowest two 
densities, our Dingle plots are ill-defined. From numerical 
simulations we deduce that even there a T-dependence of 
Tq on the scale of T/Tp will at most generate a 3% error 
in m*. Therefore, the effect of a T-dependent Tq will have 
an insignificant impact on m* in fig. 2(a). 

In conclusion, we have performed high precision mea- 
surements of the electronic effective mass in an ultra-high 
quality, tunable, two-dimensional electron system over a 
wide rage of r^, 6 > > 0.8. Performing various cross 
checks and exploring several sources for error we con- 
vinced ourselves that our data provide an accurate mea- 
surement of the density dependent impact of electron- 
electron interactions on the electron mass in a 2D system. 
Over wide stretches of density this mass rcnormalization 
can be negative. Various theoretical calculations repro- 
duce sections of our data quite well but none shows good 
agreement for the whole range of densities. 
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